Network Neuroscience

Brain network analysis has emerged as a critical framework for understanding the complex organization and function of the human brain, underpinning insights into cognition, behavior, and neuropsychiatric conditions. Central to this approach is the parcellation of the brain into discrete regions, which simplifies high-dimensional connectome data and facilitates the investigation of network architectures. However, the proliferation of brain parcellation schemes introduces significant challenges: different parcellations often yield varying network sizes and measures, complicating cross-study comparisons and the reproducibility of findings. Moreover, most connectome construction pipelines are rigid, typically outputting connectivity matrices from only one or a few parcellation schemes, which limits flexibility. In this paper, we address these issues by introducing BridgeBP, a novel toolbox designed to bridge brain parcellations by leveraging continuous brain connectivity concepts. BridgeBP transforms structural connectivity matrices derived from one parcellation scheme into matrices corresponding to more than 40 alternative schemes, standardizing analyses and enhancing the robustness of network studies. Through extensive evaluations, we demonstrate that BridgeBP enables consistent network comparisons across diverse parcellation frameworks, paving the way for more reproducible and generalizable insights in brain connectome research.

Network analysis of human brain connectivity provides a fundamental framework for identifying the neurobiological mechanisms that cause cognitive variations and neurological disorders. However, existing diagnostic models often treat structural connectivity (SC) as a fixed or optimal topological scaffold for functional connectivity (FC). This consequently overlooks the higher-order dependencies between brain regions that are critical for characterizing pathological alterations. Moreover, the distinct spatial organizations of SC and FC complicate their direct integration, as naive alignment methods may distort the inherent nonlinear patterns of brain connectivity. To address these limitations, we propose the Graph Diffusion Optimal Transport Network (GDOT-Net), which models disease-related topological evolution and achieves precise alignment between SC and FC. Unlike existing diffusion studies, the proposed model introduces an evolvable brain connectome modeling approach to infer the complex topological structure of brain networks, unveiling higher-order connectivity patterns linked to specific neuropsychiatric disorders. Furthermore, GDOT-Net incorporates a Pattern-Specific Alignment mechanism, leveraging optimal transport to align structural and functional topological representations in a geometry-aware manner. To capture nonlinear topological relationships between brain regions, a Neural Graph Aggregator Module was developed, which adaptively learns complex node interaction patterns in brain networks. By leveraging this module, GDOT-Net generates highly discriminative representations that form a robust basis for the precision diagnosis of brain disorders. Experiments on REST-meta-MDD and ADNI demonstrate that GDOT-Net surpasses SOTA methods in uncovering structural-functional misalignments and disorder-specific subnetworks. The source code is publicly available at this Link.

Linking structural connectivity (SC) to functional connectivity (FC) through mechanistic models remains challenging in network neuroscience. In this study, empirical data of diffusion magnetic resonance imaging (dMRI) and resting-state functional MRI (rs-fMRI) were used to reconstruct SC and FC of a swine connectome. We evaluated a structurally constrained Kuramoto phase-oscillator framework to reproduce resting-state FC and then assessed the models sensitivity to traumatic brain injury (TBI) and its longitudinal progression post-TBI. A joint tuning procedure was implemented to calibrate data-informed natural frequencies and global coupling strength. The tuned Kuramoto model was then used to evolve oscillator phases constrained by the SC, followed by a Balloon-Windkessel hemodynamic model. The optimized model produced significant edge-wise correspondence between averaged simulated FC and the empirical FC (r = 0.61, p < 0.001). Graph-theoretical analysis across network densities (30-50%) showed strong agreement for global efficiency, characteristic path length, and clustering coefficient, while modularity and small-worldness exhibited deviations. Longitudinal analysis of the swine TBI dataset revealed modest reductions in structure-function coupling over time but no significant differences across injury severities. These results demonstrate that optimized Kuramoto models can reproduce key functional network features while preserving inter-subject variability.

Functional connectivity in network neuroscience is traditionally characterized using time-averaged correlations between brain regions. While these summaries capture stable large-scale organization, they do not fully reflect the temporal structure of moment-to-moment interactions. Here, we investigate how the order of interaction used to represent brain dynamics shapes the organization recovered from neural data. We compare three interaction representations of fMRI dynamics: regional activation (node time series), pairwise co-fluctuations (edge time series), and higher-order triplet interactions (triangle time series); within a common topological framework using Mapper from topological data analysis (TDA). Across task and resting-state data, Mapper representations derived from pairwise co-fluctuations more distinctly segregate task conditions than activation-based or higher-order representations. This organization reflects structured coordination patterns beyond activation polarity and is driven by high-amplitude interaction events. Beyond task states, modularity quality computed across all Mapper representations is highest for edge time series and selectively associated with stable individual differences: higher modularity relates to higher conscientiousness and lower internalizing and externalizing symptom dimensions. Together, these findings suggest that behaviorally relevant information is reflected in the topology of moment-to-moment brain interactions. Topological analysis of interaction-level dynamics therefore provides a complementary and interpretable framework for linking large-scale neural coordination to cognition, personality, and mental health.

A fundamental goal in neuroscience is to understand how the brains physical architecture supports complex functional dynamics. While the relationship between structural connectivity and pairwise functional connectivity has been extensively studied, the anatomical basis of higher-order interactions remains poorly understood. In this study, we use multivariate information theory -specifically the O-information- to investigate how the human connectome constrains subsets of brain regions characterized by predominantly redundant or synergistic information sharing. By analyzing the topology and community embedding of these subsets, we reveal two different structural profiles. Redundant subsets are characterized by high internal connection density and strong weights. Their nodes have high clustering and occupy globally less central positions. In contrast, synergistic subsets consist of globally central nodes with high betweenness centrality. We further demonstrate that leveraging these structural features, in particular node centrality, significantly improves the identification of synergistic subsets compared to random sampling. Together, these results demonstrate that the human connectome imposes specific constraints on higher-order information sharing, extending structure-function relationships beyond pairwise interactions and providing new insight into the structural origins of multivariate functional organization.

The connectivity of brain networks is extremely sparse due to metabolic, physical and spatial constraints. Although wiring sparsity can confer computational advantages for biological and artificial neural networks, sparse networks require fine parameter tuning and exhibit strong sensitivity to perturbations. How brains achieve their efficiency and robustness is unclear. Here we addressed this by analysing the dynamical properties of Echo State Networks with wiring based on the Drosophila melanogaster fruit fly connectome, compared with sparsity-matched random-wiring networks. We evaluated these networks on a set of eight cognitive tasks, and found that connectome-based neural networks (CoNNs) typically showed narrowly distributed task engagement across their neurons. The importance of a neuron for task performance correlated with its node degree, local clustering, and selfrecurrency, and these correlations were stronger in CoNNs than in random networks. CoNNs were more robust to neuronal loss, retaining their task performance and beneficial dynamical properties such as criticality and spectral radius better than random networks. Similarly, CoNNs were more robust to hyperparameter variations in both input and recurrent weight scaling. Using theoretical arguments and numerical simulations, we show that excess CoNN node self-recurrency is sufficient to explain this enhanced robustness. Overall, these results identify non-random features of connectome wiring that allow brains to reconcile extreme sparsity with reliable computation. SignificanceBrain networks support robust computation even though they operate under extreme wiring sparsity due to metabolic and spatial constraints. While sparse networks typically require fine-tuning and are sensitive to perturbations, we show that biological connectomes support specialised, efficient task engagement and remain robust to neuron loss and parameter variation. We identify excess neuronal selfrecurrency as a key structural feature underlying this stability. These results reveal how non-random connectivity stabilises computation in extremely sparse networks, providing principles for understanding brain function and designing robust, efficient artificial neural systems.

Across the last three decades, functional magnetic resonance imaging (fMRI) research - through both resting-state (rsfMRI) and task-based (tfMRI) studies - has greatly advanced our understanding regarding the neural basis of cognition. Yet the mechanistic relationship between rsfMRI and tfMRI is still poorly understood. In particular, it remains unclear how and why the brain activation patterns observed during the resting state are linked to cognitive functioning and individual differences present during task performance. Here, we test a unifying computational account which postulates that task contexts modulate the nonlinear attractor landscape and associated dynamical properties of the brain present under resting conditions, and further that the nature of this modulation is impacted by meaningful cognitive individual differences. To test this account, we develop a joint rsfMRI-tfMRI modeling and analysis framework called Mesoscale Individualized NeuroDynamics with eXogenous inputs (MINDy-X) and apply it to resting and N-back working memory task data from the Human Connectome Project. We first validated that the joint model can simulate and predict both rsfMRI and tfMRI data accurately, consistent with a common underlying dynamical system. Analyses of this joint model revealed that task-related modulation bifurcated the predominantly multistable attractor dynamics present during the resting state towards a predominantly monostable dynamics observed during N-back task states. This topological shift was also accompanied by a geometric reconfiguration, with the task state characterized by an enrichment of dynamical attractor "motifs" clustered around the frontoparietal (FPN) and default mode (DMN) networks. Task-related modulations of this attractor landscape were further subject to clear individual differences, such that individuals who did not exhibit a shift in attractor topology were more error-prone and less cautious in responding, while closer geometric proximity to the FPN and DMN motifs explained additional aspects of task performance. N-back behavior was best characterized by the combination of topological and geometric properties present in both task and rest states, suggesting that they each account for unique aspects of individual variability. The current work supports a novel computational framework for understanding the whole-brain neural activity patterns observed during rsfMRI and tfMRI as reflecting different states within a common non-linear dynamical system. This framework provides a new vocabulary for characterizing cognitive functioning in terms of the unique geometric and topological configuration of the associated attractor landscapes, with the potential for wide application in many domains of basic and clinical neuroscience research.

Functional brain networks are conventionally constructed using measures of direct temporal synchrony between neural signals, implicitly restricting connectivity to scale-specific interactions. Here, we introduce an alternative framework in which interregional similarity is defined through correlations between multiscale entropy (MSE) profiles, enabling network construction based on scale-dependent dynamical structure rather than instantaneous alignment. Using resting-state fMRI data from the Human Connectome Project (N = 1003), we systematically compare MSE-based networks with conventional time-series-based networks across conventional/spectral graph-theoretical, and information-theoretic measures. We show that MSE-based networks exhibit stronger modular organization, enhanced local segregation, and distinct global integration patterns, reflecting a reorganization of functional architecture when multiscale dynamics are taken into account. Importantly, MSE-based networks demonstrate substantially greater sensitivity to biologically meaningful variability, revealing robust and reproducible sex differences across multiple network measures, in contrast to the limited and inconsistent effects observed in conventional networks. These findings suggest that multiscale representations provide a more informative and biologically grounded basis for functional brain network construction, capturing aspects of neural organization that are not accessible through direct synchrony alone.

White matter structural connectivity constrains large-scale brain communication, yet most network models do not account for biologically meaningful differences between connections. Although axonal diameter and myelination influence neural signaling at the microscale, how these features shape systems-level functional connectivity remains unclear. Here, we test whether structural connectomes weighted by white matter microstructure give rise to distinct communication regimes that differentially predict multimodal functional connectivity. Combining quantitative MRI and advanced diffusion modeling, we constructed whole-brain networks weighted by tract caliber and multiple myelin-sensitive measures. To these, we applied routing- and diffusion-based communication models and used the resulting communication metrics to predict haemodynamic and frequency-resolved electromagnetic connectivity. Myelin-weighted networks preferentially enhanced long-range communication efficiency and redistributed spectral energy toward globally integrative topological eigenmodes. In contrast, caliber-weighted networks emphasized mesoscale organization and short-range communication. Across nested regression models controlling for geometric embedding and network topology, myelin-sensitive communication explained unique variance in functional connectivity with effects varying systematically across cortical systems and frequency bands. The strongest coupling was observed for alpha-band connectivity in association and attentional networks, consistent with a role for myelin-dependent communication delays in supporting long-range alpha synchrony. These findings demonstrate how distinct white matter microstructural features give rise to heterogeneous large-scale communication regimes: tract caliber and myelin bias communication toward locally specialized and globally integrative architectures, respectively. By integrating biologically informed connectomics with communication modeling and multimodal functional data, this work advances a mechanistic account of how white matter microstructure shapes macroscale brain dynamics.

Brain decoding from fMRI data using artificial neural networks traditionally operates at the regional level, identifying which brain areas activate during tasks but ignoring how these regions interact through structural networks. While Graph Neural Networks can capture connectivity, they require prohibitively large datasets for typical neuroscience studies. We introduce a message-passing mechanism that allows a shallow neural network to incorporate structural connectivity, enabling network-level interpretation from limited data. Using motor task data from 30 Human Connectome Project subjects, we evaluate seven structural connectivity matrices derived from deterministic and probabilistic tractography. Our approach achieves 83.0% classification accuracy while revealing functional network organization. We demonstrate that sparser, anatomy-driven connectivity matrices outperform dense alternatives, and that normalizing for network size improves model performance. Critically, our method is capable of exposing structural pathways contributing towards classification, distinguishing between complete network recruitment and selective regional activation. This approach bridges the gap between high-performance brain decoding and biological fidelity of the model, enhancing neuroscientific understanding, with implications for analyzing network dysfunctions in neurological disorders such as Alzheimers disease (AD), attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder, mild cognitive impairment (MCI), and schizophrenia.

The routing of information across the brains structural network is central to its wide range of functional capabilities. However, the mechanisms underlying information routing in complex brain networks, particularly between regions that do not share a direct anatomical connection, remain poorly understood. Neural mass models (NMMs), a computational modelling framework capable of capturing complex neural dynamics across scales, can potentially be used to study the dynamical and network bases of these vital polysynaptic routing processes. In this study, we investigate polysynaptic signalling in three widely used NMMs, obeying Ornstein-Uhlenbeck, Stuart-Landau, and Jansen-Rit dynamics, by tracking the propagation of a discrete, focal, high-amplitude perturbation across the underlying network. We find that polysynaptic propagation emerges in all tested NMMs when configured within dynamical regimes that effectively enhance the persistence of perturbations. We also find distinct parameter domains that maximise signal propagation to directly connected regions and to those separated from the source by at least two hops. Finally, we benchmark in silico stimulus propagation in the brain network against an empirical dataset of direct electrical stimulation trials, to explore the relative capabilities of the NMMs in capturing signal propagation to connected versus unconnected regions. This analysis highlights the significance of dynamical repertoire in capturing stimulus propagation outcomes. Overall, this study provides insights into how dynamical and network features shape signal propagation over complex brain networks.

Resting-state functional connectivity (RSFC) and networks (RSNs) provide insight into large-scale brain organization and its disruption in neurological disease. RSNs are most commonly assessed using fMRI, yet its translational use is constrained by high cost, motion sensitivity, and limited feasibility for repeated measurements. Functional near-infrared spectroscopy (fNIRS) offers a portable alternative, but its reliability for RSFC and RSN mapping remains insufficiently established. Near whole-head fNIRS data and fMRI-BOLD signals of corresponding cortical regions were extracted, based on which RSN organization was compared across two independent cohorts of 31 participants each. Cross-modal convergence and divergence were assessed using bivariate and partial correlations across multiple network levels. Edgewise analyses revealed substantial modality differences with bivariate correlations (50-61% of edges), which were markedly reduced using partial correlations (<3%). Group-level connectivity patterns showed moderate cross-modal similarity (r {approx} 0.37). At nodal level, net strength, local efficiency, and path-length differed substantially between modalities, while normalized strength and assortativity were largely comparable. Across nodes, group-level graph-metric distributions were broadly similar for normalized strength, assortativity, local efficiency, and path length (rho {approx} 0.27-0.5). At network-level, fNIRS-derived modules significantly overlapped with fMRI modules, particularly based on bivariate correlations, identifying default mode, attentional, executive, salience, sensorimotor, and visual networks (Jaccard {approx} 0.27-0.5). Overall, fNIRS captured key features of large-scale RSFC and RSN organization observed with fMRI, supporting meaningful cross-modal correspondence and translational utility. While partial correlations enhanced edge-level agreement, they attenuated nodal and modular recovery, suggesting greater suitability for targeted connectivity analyses rather than whole-network characterization.

Analysing higher-order interactions among simultaneously recorded neurons can provide crucial insights into neural network dynamics. Recent technological advances have enabled large-scale, long-term neuronal recordings, but analysis of such datasets often relies on simpler statistics due to computational and statistical challenges in assessing higher-order interactions. Here, we developed CHOIR, an efficient and reliable method for calculating higher-order interactions from large-scale neuronal recordings. We then used the inferred HOIs to uncover the underlying functional connectivity, differentiating between connectivity motifs in the space of pairwise and triplet-wise interactions. We found that this approach could successfully distinguish stationary and running states, sleep and awake states, and neuronal ensembles with distinct activity patterns in mice. Furthermore, we identified potential circuit architectures underlying different higher-order interactions, which we confirmed through simulations of large-scale spiking networks with specific subnetwork connectivity. Applying CHOIR to a causal manipulation dataset further confirmed the role of lateral inhibition, a key inhibitory motif, in generating specific HOI patterns. Our work provides a systematic analysis of higher-order interactions in diverse datasets and suggests that HOIs can reveal circuit motifs underlying neural dynamics across brain areas and brain states.

Neurodegenerative diseases such as Mild Cognitive Impairment (MCI), Multiple Sclerosis (MS), Parkinsons Disease (PD), and Amyotrophic Lateral Sclerosis (ALS) are becoming more prevalent. Each of these diseases, despite its specific pathophysiological mechanisms, leads to widespread reorganization of brain activity. However, the corresponding neurophysiological signatures of these changes have been elusive. As a consequence, to date, it is not possible to effectively distinguish these diseases from neurophysiological data alone. This work uses Magnetoencephalography (MEG) resting-state data, combined with interpretable machine learning techniques, to support differential diagnosis. We expand on previous work and design a Riemannian geometry-based classification pipeline. The pipeline is fed with typical connectivity metrics, such as covariance or correlation matrices. To maintain interpretability while reducing feature dimensionality, we introduce a classifier-independent feature selection procedure that uses effect-sizes derived from the Kruskal-Wallis test. The ensemble classification pipeline, called REDDI, achieved a mean balanced accuracy of 0.81 ({+/-}0.04) across five folds, representing a 13% improvement over the state-of-the-art, while remaining clinically transparent. As such, our approach achieves reliable, interpretable, data-driven, operator-independent decision-support tools in Neurology.

Resting-state functional MRI has been widely used to study brain connectivity, yet the test-retest reliability of commonly used metrics remains a concern. To improve reliability, extended scan lengths and larger subject cohorts are often recommended. However, these solutions can be impractical and pose challenges, particularly in studies of clinical populations. Here, we systematically assess the reliability of two main types of functional connectivity measures: node-based connectome metrics (edge-level intraclass correlation coefficient [ICC], connectome-level ICC, functional connectivity fingerprinting, and discriminability); and voxel-based resting-state networks (RSNs) (spatial similarity of independent component analysis [ICA]-derived RSN maps quantified using the Dice coefficient). Using data from the Human Connectome Project, we evaluated the effects of scan length (3.6, 7.2, 10.8, and 14.4 minutes) and number of participants (n = 10, 20, 50, and 100), on both within-session and between-session reliability. We found that multivariate connectome metrics demonstrated greater reliability than edge-level measures, and that scan length had stronger influence on test-retest reliability than the number of participants. For connectome metrics, 14 minutes of scanning and a cohort of approximately 20 participants were sufficient to achieve reliable estimates. In contrast, RSN measures benefited from larger cohort sizes. Our findings provide practical guidelines for designing resting-state fMRI studies in terms of scan length and number of participants, balancing reliability and feasibility. Ultimately, protocol choices should be guided by the specific study objectives and the functional connectivity metric of interest.

Functional magnetic resonance imaging (fMRI) data are inherently complex, characterized by high dimensionality, intricate inter-regional dependencies, and substantial individual variability across experimental paradigms. Traditional linear mixed models (LMMs) provide a principled framework that models population-level fixed effects while estimating variance components arising from subject-level random effects; however, they often fail to adequately capture nonlinear relationships inherent in neuroimaging data. To address these limitations, we introduce the nonlinear mixed model (NMM) approach, an innovative extension of the LMM framework that integrates neural networks to flexibly model complex fixed-effect relationships while preserving the random-effects structure to account for individual differences. NMM advances fMRI analysis by: (1) identifying robust functional connectivity (FC) patterns consistently observed across multiple paradigms; (2) leveraging SHapley Additive exPlanations (SHAP) analysis to provide post-hoc interpretability of the nonlinear fixed effects, quantifying how age, sex, and paradigm contribute to predicted FC and how these effects are distributed across large-scale brain networks; and (3) using subject-specific random effects as neural fingerprints that not only show systematic variability across attention and default mode systems but also predict standardized cognitive scores, demonstrating biological relevance. Applied to the Philadelphia Neurodevelopmental Cohort (PNC) across emotion, n-back, and resting-state paradigms, NMM achieved superior model fit relative to classical LMMs, as evidenced by lower mean squared error (MSE) in predicting FC. This framework offers a statistically rigorous and practically explainable approach for modeling large-scale FC from modest covariates while explicitly separating population-level effects from stable individual variability in functional brain organization.

Biological brain networks flexibly reconfigure functional connectivity between integration and segregation through neuromodulatory systems--noradrenaline (NA) and acetylcholine (ACh)--without altering structural connectivity. Inspired by this mechanism, we propose a modular echo state network (ESN) with context-dependent NA and ACh gain modulation, where NA promotes inter-module integration via response gain and ACh promotes intra-module segregation via multiplicative gain. We evaluate the model on two context-dependent tasks: a segregation/integration task and a context-dependent decision task. Across both tasks, the modulated model consistently outperformed the baseline, with task-appropriate modulation profiles emerging naturally from optimization and functional connectivity analysis confirming context-appropriate dynamic reorganization. These results demonstrate that neuromodulatory gain control enables adaptive, context-sensitive computation in structurally fixed reservoir networks.

Neural information flow describes the movement of activity between neurons or brain areas. Advances in experimental methods have allowed production of large amounts of observational data related to neuronal activity from the single-neuron to population level. Most current methods for analysing these data are based on pairwise comparison of activity, and fall short of reliably extracting neural information flow network structure. Dynamic Bayesian networks may overcome some of these limitations. Here we evaluate the performance of a range of Bayesian network scoring metrics against the performance of multivariate Granger causality and LASSO regression for their ability to learn the connectivity underlying simulated single-neuron and neuronal population data. We find that discrete dynamic Bayesian networks are the best performing method for single-neuron data, and perform consistently for neural-population data. Continuous dynamic Bayesian networks have a tenancy to learn overly dense structures for both data types, but may have utility in scoping studies on single-neuron data. Multivariate Granger causality is the most robust method for learning structure of neural information flow between neural-populations, but performs poorly on single-neuron data. Significance testing within multivariate Granger causality produces variable results between data types. Overall, this work highlights how the analysis of neural information flow can vary depending on they type and structure of underlying data, and promotes discrete dynamic Bayesian networks as a useful and consistent tool for neural information flow analysis.

The primary objective of system neuroscience is to understand the functional mapping and its causation in the dynamics of the brain network. Some experimental and methodological studies suggest that functional modularity and its hierarchical information processing in the brain network are crucial to understanding the functional role of task-specific or state-specific information flow in the brain. However, because most of the established techniques for detecting effective network structures in the neuroscience research field are strongly based on the "Granger causality" perspective, existing causal discovery methods specified for brain network analysis cannot identify the causal hierarchy in the modular network in the brain due to spurious correlation issues and indistinguishability of causal direction under the Gaussianity of observational noise in a linear system. To address the issues, we developed a causal discovery method for synchronous neural dynamics, called the Jacobian-informed linear non-Gaussian acyclic model, "j-VAR-LiNGAM", by incorporating the information of the Jacobian matrix determined from a phase-coupled oscillator model estimated from observed neural data into the VAR-LiNGAM algorithms. The method was validated by showing that it could extract causal ordering in both synthetic data and empirical neural observed data. Moreover, by analyzing the observed neural oscillatory signals obtained from mice and humans, we confirmed that our method identified causally hierarchical structures in the brain, which aligned with the neurophysiological interpretations. These findings suggested that our proposed method can reveal the neural basis of hierarchical information processing in the brain network.

Conventional connectome edge weights, such as number of streamlines (NOS) or diffusion tensor imaging (DTI) metrics, lack specificity to microstructural details which may hold relevance for macroscale brain organisation. Since biophysical diffusion modelling offers greater specificity to microstructure, we investigated whether parameters from the Standard Model of diffusion in white matter provide informative alternatives for connectome weights - namely the intra-axonal signal fraction (f) and perpendicular extra-axonal diffusivity [Formula], as proxies of axonal density and myelination, respectively. Using diffusion MRI data from healthy adults, we constructed structural networks at four parcellation scales, weighted by f, [Formula], NOS, fractional anisotropy (FA) and radial diffusivity (RD). While all weights reproduced expected small-world properties, only [Formula] and normalised NOS captured non-random properties of local organisation across all spatial scales. We then correlated each weighted connectome with resting-state fMRI functional connectivity and intracranial measurements of conduction velocity. At the whole-brain level, although NOS gave strongest coupling with fMRI functional connectivity, only [Formula] exhibited significant structure-function coupling across all spatial scales and modalities. At the regional level, [Formula] and RD gave highest consistency in structure-function coupling across spatial scales. Thus, connectome weights derived from [Formula] capture meaningful aspects of brain network organisation with functional relevance.